AP -  A  1 4  7  485 
UNCI  ASS  1  F  I  Ell 


NONLINEAR  WAVES 'IN  A  ROD  RESULTS  FUR  INCOMPRESSIBLE 
EIASIIC  MATER  I ALS I U )  BALLISTIC  RESEARCH  LABS  ABERDEEN 
PROVING  GROUND  ML)  T  W  WRIGHT  OCT  R4  BRL-MR-3390 
SR  I  At) -F  300  491  F/G  20/11 


m 


NL 


1.4 


MICROCOPY  RESOLUTION  TEST  CHART 
NATIONAL  BUREAU  OF  STANDARDS- 1963- A 


f  SO  CD 


,  e3ooy?/ 


MEMORANDUM  REPORT  BRL-MR-3390 


NONLINEAR  WAVES  IN  A  ROD:  RESULTS  FOR 
IMCOMPRESSIBLE  ELASTIC  MATERIALS 


Thomas  W.  Wright 


October  1984 


DTIC 

ELECTEE! 
NOV  1  5  1984 

i 

B 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED. 


US  ARMY  BALLISTIC  RESEARCH  LABORATORY 

ABERDEEN  PROVING  GROUND,  MARYLAND 


Destroy  this  report  when  it  is  no  longer  needed. 
Do  not  return  it  to  the  originator. 


Additional  copies  of  this  report  may  be  obtained 
from  the  National  Technical  Information  Service, 

U.  S.  Department  of  Commerce,  Springfield,  Virginia 


22161. 


The  findings  in  this  report  are  not  to  be  construed  as  an  official 
Department  of  the  Army  position,  unless  so  designated  by  other 
authorized  documents. 


The  use  of  trade  names  ot  manufacturers'  names  in  this  report 
does  not  constitute  indorsement  of  any  commercial  product. 


SECURITY  CLASSIFICATION  of  This  PAGE  (When  Dmtm  Entmrmd) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

1.  REPORT  NUMBER  2.  GOVT  ACCESSION  NO. 

MEMORANDUM  REPORT  BRL-MR-3390  A ^7 

3.  RECIPIENT'S  CATALOG  NUMBER 

4.  TITLE  (end  Submit) 

Nonlinear  Waves  in  a  Rod:  Results  for 
Incompressible  Elastic  Materials 

5.  TYPE  OF  REPORT  8  PERIOO  COVERED 

Memorandum  Report 

6.  PERFORMING  ORG.  REPORT  NUMBER 

7.  AUTHORO) 

Thomas  W.  Wright 

8.  CONTRACT  OR  GRANT  NUMBER^ 

9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

US  Army  Ballistic  Research  Laboratory 

ATTN:  AMXBR-TBD 

Aberdeen  Proving  Ground,  MD  21005-5066 

10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  8  WORK  UNIT  NUMBERS 

1L161101AH43 

It.  CONTROLLING  OFFICE  NAME  ANO  ADDRESS 

US  Army  Ballistic  Research  Laboratory 

ATTN:  AMXBR-OD-ST 

Aberdeen  Proving  Ground,  MD  21005-5066 

12.  REPORT  DATE 

October  1984 

13.  NUMBER  OF  PAGES 

31 

14.  MONITORING  AGENCY  NAME  A  ADDRESSf/f  dlllmrmit  (torn  Controlling  Ollleo) 

IS.  SECURITY  CLASS,  (o t  thlm  report) 

IS*  DECLASSIFICATION/ DOWN  GRADING 
SCHEDULE 

IS.  DISTRIBUTION  STATEMENT  (ol  thlo  Roport) 

Approved  for  public  release;  distribution  is  unlimited. 

17.  DISTRIBUTION  STATEMENT  (ol  tho  abolroct  ontorod  In  Block  20,  It  dltloronl  from  Roport) 

18.  SUPPLEMENTARY  NOTES 

19.  KEY  WORDS  (Continue  on  rtrorto  o  Ido  If  nocoooory  end  Idonllty  by  block  number) 

steady  waves,  rod,  shock  waves,  incompressible  material,  nonlinear  elastic 

20.  ABSTRACT  CCoBthruo  on  rororoo  ohbr  ft  nocoooory  mod  Idonllty  by  block  number) 

A  nonlinear  intrinsic  theory  is  used  to  describe  the  motions  of  a  straight 
round  elastic  rod  including  the  influence  of  radial  shear  and  inertia.  Con¬ 
sideration  of  steady  wave  motions  reduces  the  two  coupled  partial  differential 
equations  to  ordinary  differential  equations  for  which  two  integrals  of  the 
motion  may  be  found.  For  incompressible  elastic  materials  with  the  restric¬ 
tion  of  small  strain  gradients,  but  arbitrary  finite  strains,  a  large  variety 
of  exact  solutions  may  be  found  by  quadrature.  These  include  large  amplitude 

DO  1  jSJTti  1473  «»T»o»  OF  *  mov  as  is  obsolete  Unclassified  (OVER) 


SECURITY  CLASSIFICATION  OF  THIS  RASE  ( When  Dolt  Enlorod) 


periodic  waves  (which  may  contain  shocks),  solitary  waves,  and  in  some  cases 
waves  that  are  transitional  from  one  stress  level  to  another.  Such  solutions 
may  be  found  for  uniform  stress/strain  curves  that  are  concave  up  or  down,  that 
contain  inflections,  and  even  curves  that  represent  phase  transitions. 


Unclassified 


MCUftITV  CLASSIFICATION  OF  THIS  FAOCfWhM  OM  Enfnd) 


TABLE  OF  CONTENTS 


LIST  OF  ILLUSTRATIONS .  5 

I.  INTRODUCTION  .  7 

II.  INCOMPRESSIBLE  MATERIAL .  8 

III.  STEADY  WAVES  .  8 

IV.  SMALL  SURFACE  ANGLE .  9 

V.  SHOCK  WAVES . 10 

VI.  STEADY  WAVE  SOLUTIONS . 15 

A.  T(e)  is  concave  down  and  at  e  ,  p(e  )  <  T  (e  ) . 15 

O  O  GO 

B.  T(e)  is  concave  down  and  at  eo,  T^e^)  <  pCe^) . 18 

C.  T(e)  has  an  inflection  point  and  TE(e0)  <  p(eo) . 18 

VII.  UNBOUNDED  SOLUTIONS . 21 

REFERENCES . 23 

DISTRIBUTION  LIST . 25 


DT1C 


B 


Accession  For  y  | 

NT IS  GRA&I 
DTIC  T.'B 

Itumnoiino  *r*3 

Jus  t  i  !  .  c: :: 

By  -  . . . 

Dint  r.-<*  A--;:/ 

— 

Co.  1;.3 

Aj  ill  •  • lit- 

;.»v 

Dial  ■  :  ,  •  , .  . 

.  i| 

LIST  OF  ILLUSTRATIONS 


FIG.  NO. 

0 

1  Stress-strain  curve  and  the  equal  area  construction 

2 

for  bounded  smooth  waves  where  vi(eQ)  <  pc  <  Te(G0)  •  • 

2  The  shear  modulus  y(e)  is  assumed  to  be  either  mono- 
tonically  decreasing  or  monotonically  increasing  .  .  . 

3  Stress-strain  curve  and  the  geometry  of  shock  waves. 

See  text,  Section  5  . 

4  Typical  mode  shapes  of  steady  waves.  Dashed  lines 
indicate  shock  waves.  Roman  numerals  refer  to 

cases  discussed  in  Section  6  . 

5  Stress  strain  curve  and  the  equal  area  construction 

2 

for  the  case  T£(eo)  <  pc  <  u(eQ)  .  Typical  mode 

shapes  are  also  shown.  Dashed  lines  indicate  shock 
waves .  . 

6  Stress  strain  curves  with  inflections  and  the  equal 

area  construction.  Case  a)  shows  a  simple  inflection; 
Case  b)  shows  a  model  for  a  phase  change.  Typical  mode 
shapes  are  also  shown  . 

7  Phase  plane  showing  both  bounded  and  unbounded  trajec¬ 

tories.  Typical  mode  shapes  for  some  of  the  unbounded 
trajectories  are  also  shown  . 


Page 

11 

13 

14 


16 


19 


20 


22 


5 


I.  INTRODUCTION 


A  rod  may  Ije  represented  intrinsically  as  a  stretching  line  endowed 
with  structure.  In  the  most  elementary  version  of  this  idea  one  scalar 
internal  variable  can  be  used  to  simulate  the  effects  of  finite  transverse 
dimensions  in  a  straight  cylindrical  rod  that  undergoes  only  axisymmetric 
motions.  The  rod  then  has  two  kinematically  independent  degrees  of  free¬ 
dom.  In  cylindrical  coordinates  R,0,Z  the  motion  of  the  rod  is  assumed  to 
be 


z  =  Z  +  w(Z,t),  r  =  R [ 1  ♦  u(Z,t)],  9*0,  (1) 

where  the  initial  radius  of  the  rod  is  a,  and  t  is  time.  All  field  quanti¬ 
ties  depend  only  on  the  axial  coordinate  Z  and  time  t. 

In  a  previous  paper2  some  of  the  consequences  of  the  theory  were 
worked  out  for  an  elastic  material  with  a  strain  energy  W(c,u,q)  that 
depends  on  axial  strain  e  =  w^,  radial  strain  u,  and  gradient  of  radial 

strain  q  =  u^,  where  the  subscripts  denote  partial  differentiation.  In 

the  terms  of  general  continuum  mechanics  W  is  the  energy  density  of  a  one¬ 
dimensional  elastic  continuum  with  one  scalar  internal  variable,  and 
because  of  invariance  considerations,  it  is  required  to  be  an  even  function 

12  12 

in  q.2  With  a  kinetic  energy  density  K  =  +  Jp2ut  Hamilton's  Princi¬ 

ple  leads  to  the  following  Euler-Lagrange  equations, 

SZ  P1  Wtt  ’  ^Z  ‘  P  =  p2Utt  (2) 

where  S  =  W£,  P  =  Wu>  Q  =  W^.  In  Reference  2  it  was  3hown  that  the  forces 

S,  P,  Q  may  be  interpreted  as  averages  of  Piola-Kirchhoff  stresses,  taken 
over  a  cross-section  of  the  rod, 

s  «->  j  J  TZZdA,  P  I  J  (TrR  +  T6  0)dA  ,  (3) 

Q  -*-+  i  J  RTrZdA  , 

and  the  appropriate  kinetic  energy  densities  are  p,  =  p  where  p  is  the 

12  1 

density  of  the  rod  material  and  p£  =  ^ a  •  Equation  (2)j  describes  forces 


lAntman,  S.  S. ,  "The  Theory  of  Rods , "  Handbuch  der  Physik ,  Vol.  Via/ 2, 
Spring er-Ver lag.  New  York  1972. 

2Wright ,  T.  W. ,  "Nonlinear  Waves  in  Rods , "  in  Proa.  IUTAM  Sump-  on  Finite 
Elastioity,  D.  E.  Carlson  and  R.  T.  Shield ,  eds. ,  Martinus-Nijhoff  Publ. , 
The  Hague,  1981. 
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and  motions  in  the  axial  direction,  and  equation  (2)2  describes  forces  and 
motions  in  the  radial  direction. 


II.  INCOMPRESSIBLE  MATERIAL 

If  the  material  is  incompressible,  then  e  and  u  are  no  longer  kine¬ 
matically  independent  and  therefore  cannot  be  varied  independently  in 
Hamilton's  principle.  The  incompressibility  relationship  between  e  and  u, 
shown  in  equation  (4), 

(1  ♦  e)(l  ♦  u)2  =  1  ,  (4) 

must  be  used  as  a  side  condition  in  the  variational  principle,  and  the 
explicit  dependence  of  the  strain  energy  on  radial  strain  must  be  sup¬ 
pressed  , 

W(e,u,q)  =  W(e,q)  .  (5) 

Under  these  conditions  the  Euler-Lag range  equations  become 

Sz  -  [ X ( 1  +  u)2]z  =  p1wtt  ,  Qz  +  2X(1  +  u)(l  +  c)  -  P2utt»  (6) 

where  X  is  a  Lagrange  multiplier,  and  now  S  *  W  and  Q  =  W  .  Clearly  X  may 

^  4 

be  interpreted  as  a  superimposed  hydrostatic  (Cauchy)  pressure.  In  the 
formulation  leading  to  either  (2)  or  (6)  more  or  less  equal  weight  is  given 
to  both  radial  shear,  which  enters  through  the  dependence  of  W  on  q,  and  to 
radial  inertia,  which  enters  from  the  finite  value  for  P2. 


III.  STEADY  WAVES 

A  steady  wave  is  a  disturbance  that  travels  down  the  rod  at  a  constant 
speed  without  distortion.  To  examine  equations  (6)  for  the  possibility  of 
steady  waves  it  is  supposed  that  all  field  variables  depend  only  on  the 
combination  £  «*  Z  -  ct  ,  where  c  is  an  arbitrary  constant.  The  partial 
differential  equations  then  become  ordinary  differential  equations. 

S'  -  [X(l  +  u)2]'  =  pc2e'  (7) 

Q'  ♦  2 X C 1  ♦  u)(l  +  e)  =  \  pa2c2u"  (8) 


where  the  dash  signifies  differentiation  with  respect  to  £ ,  and  Pj  and  P2  have 
been  written  in  terms  of  p.  As  in  Reference  2  there  are  two  integrals  of 
the  motion. 


(9) 

(10) 


S  -  A(1  +  e)'1  =  pc2(  1  +  e)  +  A  , 

Qu'  +  S(1  +  e)  -  X  -  W  = 

0 

19  2  19  2 

pc  [(1  ♦  t)  +  i  a'  u'  1  ♦  B  , 

where  A  and  B  are  constants  of  integration.  The  first  integral  is  straight 
forward.  The  second  integral  is  obtained  by  multiplying  (7)  by  (1  +  e), 

(8)  by  u' ,  adding  the  two  resulting  equations,  and  noting  that  (10)  is  the 
integral  of  the  summed  expression.  Equation  (9)  may  be  used  to  eliminate  A 
from  (10),  arriving  finally  at  equation  (11), 

Qu’  ♦  A(1  ♦  e)  -  W  =  |  pc2[~  a2u'2  -  (1  ♦  t)2]  +  B  .  (11) 

So  far  the  results  are  exact  given  the  initial  premises  concerning  the 
strain  energy,  the  kinetic  energy,  and  the  condition  of  incompressibility. 


IV.  SMALL  SURFACE  ANGLF 

Now,  suppose  that  tht  surface  >"el  -  is  small  in  the  sense  that 
au'  <<  1,  and  express  the  strain  energy  as  the  first  two  terms  of  a  power 
series  in  u*  with  coefficients  that  depend  on  e  in  an  arbitrary  way, 

W(e,q)  3  W(e,0)  +  j  q2  W^(e,0)  +  -  .  (12) 


Alternatively,  consider  the  special  case  when  W  can  be  expressed  exactly  as 

W(e,q)  =  J  T(e)de  +  ^  a2p(e)q2  .  (13) 

T(e)  is  engineering  stress  (force  per  unit  original  area)  for  uniform 
extension,  and  y(e)  is  interpreted  as  the  shear  modulus.  (In  the  linear 
version  of  the  theory  p  would  be  the  bulk  shear  modulus  as  in  Reference  2, 
so  it  seems  fair  to  give  it  that  name  for  finite  extension  as  well.)  From 
the  assumed  form  of  the  energy  (13)  with  A  eliminated  by  use  of  (8)  and  u 
eliminated  by  use  of  the  incompressibility  condition  (*4),  the  integral  (9) 
may  be  written 

T  =  T(e)  -  pc2(e  -  e  )  -  yY  e’2  -  Y(Oe"  >  (14) 

O  O  £  c- 


1  2  -3  2 

where  Y  -  s-a  (1  +  e)~  [u(e)  -  pc  ] ,  and  T  =  T(e  ).  For  convenience  A 
o  2  o  u 

has  been  replaced  by  Tq  -  pc  (1  +  eQ) ,  Similarly  the  integral  (11)  may  now 

be  written 


1  2,  2,  ,2 

T  a  (u  -  pc  )u’  = 


Tde  -  [Tq  -  pc2(l  ♦  e0)J(l  ♦  e)  -  j  Pc2(l  +  e) 2  +  B. 
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By  rearranging  terms  and  using  the  incompressibility  condition  (4),  the 
integral  may  be  rewritten  as 


~  a  PC-J-  e' 2  ■  f  [T(e)  -  TQ  -  pc2(e  -  eo)]de  +  C  ,  (16) 

(1  +  e)  J eo 

where  C  is  an  arbitrary  constant.  Equation  (16)  also  follows  from  (14) 

directly  upon  noting  that  multiplication  by  e'  makes  the  terms 
13  12 

J  Yee<  *  YG* 1 £’  the  exact  derivative  of  y  y(e)e'  .  Solution  of  (16)  then 

follows  simply  by  quadrature.^  The  integral  on  the  right  hand  side  of  (16) 

is  the  area  between  the  homogeneous  extension  curve  and  a  straight  line 

2 

through  the  point  (T  ,  eQ)  with  slope  pc  .  As  shown  in  Figure  1,  the 
cross-hatched  areas  both  to  the  right  and  to  the  left  of  (Tq,  eq)  are 
positive.  The  Lagrange  multiplier  is  given  by  the  following  expression. 


r  \  3 

W  -  PC2 

„  2 
,2  w  pc  ,  , 

.2  2 
L(1  +  rl 

(1  *  e,3J 

e  +  2 

(1  +  e) 

V.  SHOCK  WAVES 

Although  integration  of  ( 1 6 )  is  straightforward,  the  characteristics 
of  a  particular  solution  depend  on  the  parameters  in  the  problem  and  the 
nature  of  the  functions  T(e)  and  y(e).  In  addition,  it  is  also  possible 
for  discontinuities  to  occur  in  certain  limiting  cases.  Therefore,  before 
sketching  out  the  rich  variety  of  solutions  available,  it  is  necessary  to 
consider  the  propagation  of  shock  waves. 

Since  the  cross  sectional  area  of  the  rod  cannot  change  discontinu- 
ously,  the  radial  3train  must  be  continuous  according  to  (1),  and  since  the 
material  is  incompressible,  the  axial  strain  must  be  continuous  as  well. 


t Aifantis  and  Semin 5  in  discussing  phase  transitions  in  the  presence  of 
surface  tension  and  Colemarf  in  developing  a  static  theory  of  necking  and 
drawing  of  fibers  have  encountered  equations  of  the  form 

T  =  t(e)  +  6(e)  ' 2  +  y(e)e".  In  fact ,  for  static  cases  when  c  =  0,  equa¬ 
tion  (14)  coincides  exactly  with  Coleman's  result  if  his  equtlibrium  equation 
is  obtained  from  a  variational  principle  (see  equation  (3.20)  and  the  follow¬ 
ing  discussion  in  Reference  6).  The  more  general  equation  was  solved  by 
finding  an  integrating  factor  (e.g.}  see  Reference  6 ,  equation  (3.3))  to 
obtain  a  first  integral,  followed  by  a  quadrature  to  obtain  a  second 
integral. 
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Stress-strain  curve  and  the  equal  area  construction  for 
bounded  smooth  waves  where  u  ?  <  t  r£  i 


However,  both  u'  and  e'  may  have  discontinuities,  and  therefore,  by  (17)  so 
can  X.  The  appropriate  Jump  condition  comes  from  the  integral  form  of  (2) 
or  ( 6  ),3 » 4 

\  pa2V2  -  ,  (18) 

where  [*]  signifies  the  jump  in  a  quantity  across  a  shock  wave Substitu¬ 
tion  for  W  from  (13)  gives 

pV2  =  u(e)  ,  [q]  is  arbitrary,  (19) 

so  that  the  speed  V  of  a  shock  wave  of  any  amplitude  is  determined  solely 
by  the  instantaneous  shear  modulus. 


If  the  shock  is  incorporated  into  a  steady  wave,  then  V  =  c,  and  if  y(e) 
is  a  monotonic  function  either  increasing  or  decreasing  (the  only  cases 
considered  in  this  paper),  then  for  a  given  steady  wave  speed  c,  (19)  can 
be  satisfied  only  at  isolated  values  of  e,  say  e  =  e*,  unless  u  is  truly 
constant.  The  situation  is  shown  in  Figure  2.  Since  the  left  hand  side  of 

2 

(16)  is  zero  at  a  shock  wave,  no  matter  what  the  value  of  e'  ,  two  steady 

waves  with  the  same  speed  may  be  joined  together  at  the  shock,  but  of 

course  each  wave  has  distinct  values  lor  Tq,  eq,  and  the  constant  C.  These 

values  are  not  independent  in  this  case,  and  in  fact  if  one  is  known  for  a 
particular  steady  wave,  the  others  may  be  easily  computed.  This  is  best 
seen  by  referring  to  Figure  3.  For  the  case  shown  e  >  e*,  the  point  R  has 
coordinates  (TqJ  eQ),  and  C  =  -  Area  (PQR).  If  a  second  steady  wave  has 

the  same  speed  so  that  Q'R'  is  parallel  to  QR,  then  R*  has  coordinates 

(r  ,  cM,  and  C*  =  -Area  (PQ’R').  The  quantity  [e'2]  may  be  determined 

from  (14),  noting  that  y  s  0  for  a  steady  shock  wave. 


2  12'Jf^e*^ 

[To  -  pc  (e*  -  e0)]  =  -  a  - - - 

(1  e*) 


3  (e’2]  • 


(20) 


The  magnitude  of  the  left  hand  side  of  (20)  is  just  the  vertical  distance 
QQ'  in  Figure  3* 


3 Nunziato ,  J.  W.}  and  Walsh ,  E.  K. ,  " One -dimensional  Shook  Waves  in  Uniformly 
Distributed  Granular  Materials Int.  J.  Solids  and  Structures  14  (1978) 
681-689.  ~ 

4 Wright,  T.  W. ,  " Weak  Shooks  and  Steady  Waves  in  a  Nonlinear  Rod  or  Granular 
Material,"  (to  appear ,  Int.  Solids  and  Structures) . 

lzn  a  three-dimensional  context }  a  discontinuity  in  the  strain  gradient  would 
be  classified  as  an  acceleration  wave y  but  in  the  present  context ,  where  the 
discontinuity  occurs  in  derivatives  of  lower  order  than  the  highest  that 
appear  in  the  governing  differential  equation  (2)  and  (8),  it  is  legitimate 
to  use  the  word  shock. 
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Figure  3.  Stress- strain  curve  and  the  geometry  of  shock 
See  text.  Section  5. 


VI.  STEADY  WAVE  SOLUTIONS 


By  a  process  similar  to  that  used  in  References  5  and  6  it  is  possible 
to  show  tha£  bounded  smooth  solutions  in  an  infinite  rod  must  be  one  of 
three  kinds:  periodic  waves,  solitary  waves  with  either  a  single  propa¬ 
gating  bulge  or  neck,  or  a  propagating  smooth  transition  from  one  state  of 
uniform  strain  to  another.  But  in  addition,  solutions  with  discontinui¬ 
ties,  as  described  in  the  previous  section,  may  occur.  Most  of  these  will 
appear  as  limiting  cases  of  the  smooth  solutions.  At  this  point  it  is 
easiest  to  proceed  by  considering  examples  and  special  cases.  Figure  4 
shows  a  sketch  of  the  rod  configuration  for  each  of  the  cases  considered. 

A.  T(e)  is  concave  down  and  at  eo>  <  T£(eo). 

Refer  to  Figure  1  and  consider  the  following  construction.  Through 
the  point  T  ,  eq  draw  a  straight  line  with  positive  slope  pc2,  chosen  such 

2 

that  w(E0)  <  pc  <  T  (eQ)  .  From  Figure  2  find  e*  ,  which  will  lie  to  the 

left  of  eQ  if  y(e)  is  a  decreasing  function  or  to  the  right  of  eq  if  y(e) 

is  an  increasing  function.  Finally  mark  off  equal  cross  hatched  areas  to 
the  right  and  left  of  eQ  with  e*  lying  outside  the  cross  hatched  interval, 

and  set  the  constant  C  equal  to  minus  one  half  the  total  cross  hatched 

2 

area.  From  the  construction,  it  is  clear  that  e'  £  0  over  the  whole 
interval  from  emin  to  emax  ,  and  e'  =  0  only  at  the  end  points.  Further¬ 
more,  the  construction  is  such  that  near  either  end  point  e'  varies  as 
the  square  root  of  the  distance  from  the  end  point.  Every  construction  of 
this  kind  corresponds  to  a  periodic  steady  wave.  As  the  magnitude  of  C 
increases,  the  cross  hatched  areas  extend  farther  from  e  ,  corresponding  to 

a  larger  amplitude  of  the  wave,  until  one  of  the  following  finally  occurs: 

either  i)  e  =  e  ,  or  ii)  e  .  «  e*  for  the  case  of  decreasing  y  ,  or  iii) 

max  min 

emax  =  e*  ^cr  tlie  case  increasing  y  .  The  preceding  discussion  shows 

that  periodic  waves  form  a  family  with  three  parameters,  say  Tq  ,  c,  and 

e  .  T  may  be  chosen  arbitrarily,  but  c  and  e  can  only  be  chosen 
max  o  max 

within  limits. 

The  limiting  cases  require  special  treatment. 

i)  Solitary  Waves,  e  =  e.  As  e  increases  towards  e  ,  the  period 
1  ’  max  max 

of  the  wave  also  increases  and  becomes  infinite  when  e  =  e .  The  wave  in 

max 


sAif 'antis,  E.  C.,  and  S err  in,  J.  B. ,  "Towards  a  Mechanical  Theory  of  Phase 
Transformations, "  Teoh.  Rpt. ,  Corrosion  Research  Center,  University  of 
Minnesota,  Minneapolis,  1982. 

6 Coleman ,  B.  D. , "Necking  and  Drawing  in  Polymeric  Fibers  Under  Tension,"  Arch. 
Rat.  Meah.  Anal.  83  (1982)  115-137. 
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Periodic  Waves 


v  ■■■  — _ _ _ J 


i)  Solitary  Wave 

ii)  Periodic  Shocks, 

fi  Decreasing 


t 


iii)  Periodic  Shocks, 
ft.  increasing 


ii)'  Sharp  Crested 
Solitary  Wave 


iii)'  Finite  Length 

Solitary  Wave 


iva)  Similar  to  Linear 
Dispersive  Pattern 


ivb)  Structured  Shock 


Figure  4.  Typical  mode  shapes  of  steady  waves.  Dashed  lines  indicate  shock 
waves.  Roman  numerals  refer  to  cases  discussed  in  Section  6. 


this  case  consists  of  a  single  crest,  a  propagating  bulge  in  the  radius  of 
the  rod  for  the  stress/strain  curve  shown  in  Figure  1.  At  infinity  behind 
and  ahead  of  the  wave  the  strain  is  e,  decreasing  to  em^n,  as  determined  by 

the  equal  area  rule,  at  the  crest.  Since  the  maximum  strain  is  determined 
by  the  choice  of  TQ  and  c,  solitary  waves  form  a  two  parameter  family. 

ii)  Shock  Waves,  e_.  3  e*.  As  e_._  decreases  towards  e* ,  the  crests 

rain  mm 

of  the  wave  become  less  and  less  rounded  until  finally  when  e  .  =  e* ,  the 

*  min 

crests  become  sharp  with  e'  changing  discontinuously.  Since  C  in  (16)  has 
been  chosen  such  that  both  sides  vanish  when  e  =  e*. 


lim  e'  2  *» 
e  -*■  e* 


16(1  +  e*)3 

2 

a 


T(e*)  -  To  -  pc2(e*  -  eQ) 
Pe(£*) 


(21) 


and  e'  jumps  from  one  of  the  values  that  satisfies  (21)  to  the  other.  The 
numerator  in  the  second  factor  of  (21)  is  just  the  negative  of  the  vertical 
distance  PQ  in  Figure  3.  Choice  of  c  determines  e* ,  and  choice  of  Tp 

determines  emax  from  the  equal  area  rule,  so  waves  of  this  kind  form  a  two 

parameter  family. 

iii)  Shock  Waves,  emax  3  e* .  This  case  is  the  same  as  ii)  except 

that  shock  waves  occur  at  e  „„  rather  than 

max  mm 


It  may  happen  that  two  of  the  limiting  cases  occur  simultaneously. 

ii) '  Shock  Wave,  £max  3  £  and  em^n  3  e*.  In  this  case  the  solitary 

wave  has  a  sharp  crest  on  it.  The  two  cross  hatched  areas  in  Figure  1 

extend  all  the  way  from  e  to  £  and  from  to  e*.  Since  only  one  choice  of 

o  o 

T  will  make  the  areas  equal  for  a  given  c,  these  waves  form  a  one  param¬ 
eter  family. 

iii) '  e  =  e  =  e* .  By  analogy  with  the  previous  case  it  might  seem 

max 

that  this  case  would  correspond  to  a  single  shock  wave  occurring  at  the 

maximum  strain,  but  in  fact  e'  is  zero  at  This  occurs  because  the 

max  2 

right  side  of  (16)  now  has  a  double  zero  at  e  ,  but  the  coefficient  of  e' 

max 

on  the  left  side  has  only  a  simple  zero.  Thus,  these  waves  are  simply 
periodic  waves  of  maximum  wave  length.  However,  since  discontinuities  can 
occur  at  e* ,  one  cycle  of  a  wave  could  be  joined  to  a  uniform  state  both 
behind  and  ahead  of  the  wave  thus  creating  a  --'lltary  wave  of  finite 
wavelength.  Since  the  choice  of  c  determines  e*  (and  hence  and  the 

equal  area  construction  will  determine  all  other  parameters,  these  waves 
are  a  one  parameter  family. 
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iv)  As  stated  in  Section  5,  a  shock  wave  may  join  two  steady  waves 
together  if  they  have  the  same  speed.  Cases  ii),  ii)’,  and  iii)  suggest 
many  ways  in  which  this  could  be  done.  Two  sharp  crested  periodic  waves 
with  unequal  wave  lengths  could  be  joined,  for  example,  but  two  more 
interesting  cases  are  a)  a  sharp  crested  periodic  wave  Jointed  to  a  sharp 
crested  solitary  wave,  which  resembles  the  familiar  dispersive  pattern  in 
linear  elasticity  for  a  step  load  on  the  end  of  a  right  circular  cylinder, 
and  b)  a  uniform  state  joined  to  a  sharp  crested  solitary  wave,  which  is 
really  Just  a  limiting  form  of  iva)  and  might  be  termed  a  structured  shock 
wave.  With  these  constructions  it  is  possible  to  join  two  arbitrary 
(within  limits)  stress  states  with  a  steady  wave. 

B.  T(e)  is  concave  down  and  at  eQ»  Tc(eQ)  <  . 

If  the  stress/strain  curve  for  uniform  extension  turns  over  so  far 
that  T  (e  )  <  u(e  ) ,  then  the  previous  construction  for  bounded  solutions 

must  be  modified  by  choosing  c  such  that  T  (£q)  <  pc“  <  p(eq)  •  In  all 

other  respects,  however,  the  construction  is  the  same.  Equal  areas  to  the 

right  and  left  of  extend  to  £max  and  Limiting  cases  occur  when 

coincides  with  either  r.  or  e*  (for  the  case  of  increasing  P)  or  when  ^ 

coincides  with  c*  (for  the  case  of  decreasing  u).  Periodic  waves,  includ¬ 
ing  periodic  shock  waves,  look  much  the  same  a3  before,  but  the  wave  of 
infinite  period  (i.e.,  the  solitary  wave)  is  now  a  propagating  neck  rather 
than  a  bulge  as  previously  and  the  limiting  case  is  a  sharp  notch  rather 
than  a  sharp  crest.  Since  two  different  steady  waves  can  be  joined  through 
a  shock  wave,  many  combinations  are  possible,  but  perhaps  the  most  inter¬ 
esting  case  is  the  structured  shock  wave  which  joins  a  uniform  state  to 
half  of  a  sharp  notch.  The  sketches  in  Figure  5  should  be  compared  with 
those  in  Figure  4. 

C.  Tfe)  has  an  inflection  point  and  T  (e  )  <  p(e  ). 

v  '  t.  o  o 

In  this  case  the  equal  area  construction  can  lead  to  a  smooth 
transition  from  one  uniform  state  to  another  provided  that 

T  (e  )  <  pc^  <  u(e  )  and  e*  lies  outside  the  range  of  strains  in  the  tran- 

e  o  o 

sition.  The  constant  C  i3  equal  to  half  the  doubly  cross  hatched  area 
shown  in  Figure  6a.  The  construction  is  valid  even  if  the  stress/strain 
curve  for  uniform  extension  has  a  maximum  and  a  minimum  as  shown  in  Figure 
6b.  Such  a  curve  has  been  suggested  by  Ericksen  as  a  model  for  phase 
change7  and  was  used  by  Coleman  to  describe  necking  and  drawing  of  poly- 
mera6  In  the  latter  context  the  solution  represented  by  Figure  6b  shows 
the  influence  of  radial  inertia  and  radial  shear  on  the  drawing  process  at 
high  speed.  Note  also  the  implication  that  the  drawing  speed  is  limited  by 
the  shear  speed  of  the  undrawn  (drawn)  material  for  u(e)  increasing 
(decreasing) . 


7 Ericksen,  J. 


"Equilibrium  of  Bare,"  J.  Elasticity  5  (1975)  191-201. 
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Dashed  lines  indicate  shock  waves. 


It  is  possible  for  e*  to  lie  at  either  the  minimum  or  the  maximum 
strain  of  the  transition.  In  that  case  there  will  be  a  shock  discontinuity 
on  one  side  or  the  other  of  the  transition,  but  as  long  as  the  equal  area 
construction  can  be  maintained,  a  transition  solution  still  exists.  In 
this  case  the  constant  C  is  equal  to  half  the  total  cross  hatched  area  in 
Figure  6b. 


VII.  UNBOUNDED  SOLUTIONS 

So  far  only  solutions  that  are  bounded  for  -«  <  2  -  ct  <  «  have  been 
described,  but  there  are  unbounded  steady  wave  solutions  as  well.  To  find 
these,  it  is  easiest  to  adopt  a  slightly  different  point  of  view  based  on 
the  phase  plane. 

As  an  example  refer  to  Figure  1,  consider  the  case  shown  there  for  u(e) 
decreasing  and  e*  <  ep,  and  suppose  that  and  c  are  both  fixed.  The 

construction  is  such  that  the  right  hand  side  of  equation  (16)  has  a  min¬ 
imum  at  e  =  cQ  and  a  maximum  at  e  •  e,  The  constant  C  can  be  either  posi¬ 
tive  or  negative  so  the  right  hand  3ide  of  (16)  has  either  one  or  three 
real  zeros  including  a  possible  double  zero  at  the  maximum  or  minimum.  The 

2 

coefficient  of  s'  on  the  left  hand  side  of  (16)  has  only  a  simple  zero  at 
e  =  e* .  Bearing  these  facts  in  mind  and  taking  account  of  the  relative 
positions  of  the  zeros  on  both  sides,  one  may  solve  equation  (16)  for  e'  as 
a  function  of  e  and  sketch  a  family  of  curves,  as  in  Figure  7,  each  member 
of  which  corresponds  to  a  different  value  of  C.  Because  of  the  zero  at  e*  , 
no  curve  can  cross  the  vertical  line  e  *  e*  except  the  case  where  the  right 
hand  side  of  (16)  has  a  zero  at  e*  as  well,  and  because  of  incompressibil¬ 
ity,  no  curve  can  cross  the  line  e  =  -1.  For  other  values  of  Tq,  eq  ,  and  c 

the  phase  plane  diagram  can  be  substantially  modified.  For  example,  the 
location  of  e*  obviously  will  have  a  pronounced  effect,  and  if  the  straight 
line  through  Tq,  eq  has  no  other  intersection  with  T(e)  ,  there  can  be  no 

orbits  that  correspond  to  smooth  solitary  waves. 

The  closed  trajectories  around  eQ  represent  periodic  waves  with  the 

curve  labeled  L  being  the  limiting  case  of  sharp  crested  waves.  Examples 
of  unbounded  trajectories  are  labeled  A  through  G.  Even  though  these 
trajectories  are  unbounded  in  the  phase  plane,  the  corresponding  deforma¬ 
tion  patterns  are  readily  interpretable,  and  are  also  sketched  in  Figure  7, 
showing  the  rod  radius  as  a  function  of  5.  Although  no  solution  to  an 
initial/boundary  value  problem  could  include  any  point  with  e  =  -1  or  e  =  ® 
or  e'  =  ±® ,  a  segment  of  an  unbounded  trajectory  that  does  not  include  such 
singularities  would  be  perfectly  acceptable.  Problems  in  which  the  mater¬ 
ial  moves  at  constant  speed  with  respect  to  fixed  boundaries,  where  fixed 
boundary  conditions  are  specified,  would  fall  in  this  category. 
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V 


showing  both  bounded  and  unbounded  trajectories.  Typical  mode 
ome  of  the  unbounded  trajectories  are  also  shown. 
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